Mon. Not. R. Astron. Soc. 000, HHT^ (2004) Printed 2 February 2008 (MN WF^ style file v2.2) 



Stellar mass black hole binaries as ULXs 



S. A. Rappaport^, Ph. Podsiadlowski^ and E. Pfahl^ 

1 Department of Physics, MIT 37-602B, 77 Massachusetts Ave, Cambridge, MA 02139 ('sarOmit . eduj 
^ Department of Physics, Oxford University, Oxford, UK (podsiSastro . ox . ac .ukj 

^ Chandra Fellow; Harvard-Smithsonian Center for Astrophysics; 60 Garden St., Cambridge, MA 02138 ('epf aliiacfa.harvard.edu) 
2 February 2008 

ABSTRACT 

Ultraluminous X-ray sources (ULXs) with > 10'^^ ergs s~^ have been discov- 
ered in great numbers in external galaxies with ROSAT, Chandra, and XMM. The 
central question regarding this important class of sources is whether they represent an 
extension in the luminosity function of binary X-ray sources containing neutron stars 
and stellar-mass black holes (BHs), or a new class of objects, e.g., systems contain- 
ing intermediate-mass black holes (100-1000 Mq). We have carried out a theoretical 
study to test whether a large fraction of the ULXs, especially those in galaxies with 
recent star formation activity, can be explained with binary systems containing stellar- 
mass black holes. To this end, we have applied a unique set of binary evolution models 
for black-hole X-ray binaries, coupled to a binary population synthesis code, to model 
the ULXs observed in external galaxies. We find that for donor stars with initial masses 
> 10 Mq the mass transfer driven by the normal nuclear evolution of the donor star 
is sufficient to potentially power most ULXs. This is the case during core hydrogen 
burning and, to an even more pronounced degree, while the donor star ascends the 
giant branch, though the latter phases lasts only ~5% of the main sequence phase. We 
show that with only a modest violation of the Eddington limit, e.g., a factor of ~10, 
both the numbers and properties of the majority of the ULXs can be reproduced. One 
of our conclusions is that if stellar-mass black-hole binaries account for a significant 
fraction of ULXs in star-forming galaxies, then the rate of formation of such systems 
is ~ 3 X 10^^ yr^^ normalized to a core-collapse supernova rate of 0.01 yr^^. 
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1 INTRODUCTION 

The Chandra X-ray Observatory and the XMM mission have 
been used to study entire populations of accretion pow- 
ered binary X-ray sources in external galaxies. At distances 
exceeding ~10 Mpc, the sources that can be studied are 
limited to the luminous end of the distribution function, 
e.g., Lx^W'^^ ergs s^^. The observed sources with Lx up 
to a few X 10^* ergs s~^ are very likely closely related to 
the high- and low-mass X-ray binaries that have been well 
studied for the past four decades in our own Galaxy and 
its neighbors. However, the discovery of ultraluminous X- 
ray sources (ULXs) with Emstein (Fabbiano 1989), ROSAT 
(Colbert & Ptak 2002; Roberts & Warwick 2000), and ASCA 
(Makashima et al. 2000) has been greatly extended by Chan- 
dra and XMM with their far superior sensitivity (see, e.g., 
reviews by Fabbiano & White 2004; Colbert & Miller 2004). 
These sources are typically defined to have > 10^® ergs 



s~^ (2 — 10 keV) and have been observed to luminosities as 
high as a few x 10'*" ergs s~*. A key question which observa- 
tions of these sources seek to answer is whether the compact 
object is (1) a neutron star of mass ~ 1.4 Mq or black hole of 
up to ~ 15 Mq (see, e.g., Tanaka & Lewin 1995; Greiner et 
al. 2001; Lee et al. 2002; McClintock & Remillard 2004), or 
(2) a black hole of "intermediate mass", e.g., 100— 1000 Mq 
(e.g., Colbert & Mushotzky 1999). This is the question we 
address on theoretical grounds in the current work. 

ULXs appear in different types of galaxies, including el- 
lipticals (Angelini, Loewenstein, & Mushotzky 2001; Jeltema 
et al. 2003) where their luminosities are generally confined 
to < 2 X 10^^ ergs s"* (Irwin, Bregman, & Athey 2004). 
ULXs, a few with luminosities as high as ~ 5 x 10'*" ergs s~*, 
are especially prevalent in galaxies with starburst activity, 
including ones that have likely undergone a recent dynami- 
cal encounter (e.g., Fabbiano, Zezas, & Murray 2001; Wolter 
& Trinchieri 2003; Belczynski et al. 2004; Fabbiano & White 
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2004; Colbert & Miller 2004). Two highly photogenic exam- 
ples are the Antennae and Cartwheel galaxies. The Anten- 
nae galaxies (Fabbiano et al. 2001) include 49 very luminous 
X-ray sources, 18 of which are classified as ULXs (Zezas et al. 
2002), and arc likely to be by-products of the star formation 
triggered by the collision of these galEixies (e.g., Hernquist & 
Weil 1993). Approximately half of the ULXs in the Anten- 
nae are identified with young star-forming regions while the 
other half have no apparent counterpart (sec, e.g., Fabbiano 
& White 2004). The Cartwheel galaxy reveals a substantial 
number of resolved as well as some unresolved point sources 
in a ring coinciding with starburst activity and punctuated 
by numerous HII regions (Woltcr & Trinchicri 2003, 2004; 
Gao et al. 2003). This ring of star formation is apparently 
propagating outward at ~ 50 km s~^, and the original dis- 
turbance was presumably triggered by the penetration of a 
smaller galaxy some 5 x 10* years ago. The Chandra sensi- 
tivity limit at the distance of the Cartwheel (~ 120 Mpc) is 
Lj; ~ 5 X 10^** orgs s"^ 

The ULXs found in these galaxies and numerous oth- 
ers have been suggested to harbor "intermediate-mass black 
holes" (IMBHs, e.g., Colbert & Mushotzky 1999). The mo- 
tivation for this is clear. The Eddington limit for spherically 
symmetric accretion is given by 

OO M 1 1 

LEdd ^ 2.5 X 10^*^(1 + X)-^ ergs s"' , (1) 

where M is the mass of the accretor, X is the hydrogen mass 
fraction in the accreted material, and Thomson scattering is 
taken to be the dominant source of opacity. The Eddington 
limit for neutron stars is only ~ 2 x 10'** ergs s^^, though a 
few accretion-powered X-ray sources have typical persistent 
luminosities of ~ 5 — 8 x 10^* ergs s~^ (e.g., Levine et al. 
1991, 1993) — not quite in the ULX range. Intermediate-mass 
black holes, by contrast, would have Eddington luminosi- 
ties of ~ 10^'^ — 10*^ ergs s~^ which nicely cover the ULX 
range. Moreover, the expected spectra from intermediate- 
mass black holes accreting substantially below their Edding- 
ton limit (as would be the case for most of the ULXs if they 
were IMBHs), would have low inner-disk temperatures, as 
is inferred for some of the ULXs (Miller et al. 2003; Miller, 
Fabian, & Miller 2004a,b; Cropper et al. 2004). However, 
before invoking a new type of hitherto unobserved object, it 
is sensible to ascertain whether stellar-mass black holes of 
mass 5 — 15 Mq could explain many or most of the ULXs. 
The corresponding Eddington limit for these sources is as 
high as ~ 2 X 10^^ ergs s~^, i.e., extends into the low end of 
the ULX range. Several Galactic black-hole transient X-ray 
sources have been suspected of exceeding their respective 
Eddington limits, but the best case may be GRS 1915-1-105 
where the observed Lx has been above 10^^ ergs s~^ about 
~30% of the time in daily RXTE ASM averages over the 
past 8 years (A. Levine, private communication). On rare 
occasions, and for brief intervals, Lx for GRS 1915+105 can 
be high as 7 x 10'*'^ ergs s~^ (2 — 10 keV; e.g., Greiner et 
al. 1996), which is a few times greater than the Eddington 
limit for the ~ 14 Mq black hole in this system (Greiner 
et al. 2001). If taken at face value, however, sources with 
Eddington limited stellar-mass black holes would still fall 
somewhat short in explaining a significant fraction of the 
observed ULXs. 



A number of ideas have been put forth for ways to cir- 
cumvent the problem of how ~ 10 Mq black holes could 
have apparent Lx values as high as a few x 10''° ergs s~^. 
King et al. (2001) suggested that the radiation may be geo- 
metrically beamed so that the true value of Lx does not, in 
fact, exceed the Eddington limit. Kording, Falcke, & Markoff 
(2002) proposed that the apparently super-Eddington ULXs 
are actually emission from microblazer jets that are relativis- 
tically beamed along our line of sight. However, studies of 
the giant ionization nebulae surrounding a number of the 
ULXs (PakuU & Mirioni 2003) seem to confirm the full lu- 
minosity inferred from the X-ray measurements. Begelman 
(2002) and Ruszkowski & Bcgclman (2003) found that in ra- 
diation pressure dominated accretion disks super-Eddington 
accretion rates of a factor of ~10 can be achieved due to the 
existence of a photon-bubble instability in magnetically con- 
strained plasmas. They propose that this instability results 
in a large fraction of the disk volume being composed of 
tenuous plasma, while the bulk of the mass is contained 
in high-density regions. The photons then diffuse out of 
the disk mostly through the tenuous regions, thereby efi^ec- 
tively increasing the Eddington limit. This efi^ect is shown to 
grow as M^'^^ and may yield a super-Eddington factor (here- 
after "Begelman factor") of ~10 in disks around stellar-mass 
black- hole systems.^ 

In the present work we investigate whether most of the 
ULX population — at least in galaxies with current or recent 
star formation activity — is consistent with black-hole bina- 
ries of conventional mass with only mildly super-Eddington 
luminosities. In order to carry out this study we combine a 
unique grid of binary evolution models with a binary popu- 
lation synthesis code to compute theoretical X-ray luminos- 
ity functions vs. time after a burst of star forming activity 
has occurred. We then test how well these calculated results 
match the observations. To the extent that our models are 
successful, we can invert the problem and constrain some of 
the uncertain input physics in binary stellar evolution cal- 
culations. 

In §2 we describe how the incipient black-hole binaries 
are generated. We also present a grid of 52 black-hole binary 
models evolved through the phase of mass transfer onto the 
black hole, i.e., the X-ray phase. In §3 we discuss how we 
utilize the population synthesis tools to generate X-ray lu- 
minosity functions for black-hole binaries as a function of 
time since a discrete star formation event. The results are 
presented as color images of the evolving luminosity function 
vs. time, as well as quantitative line plots of the populations 
vs. X-ray luminosity and vs. time. In §4 we discuss how 
our results apply to ULXs, particularly those that are being 
observed in increasingly large numbers in external galaxies 
with Chandra and XMM. 



We note that there are essentially no empirical guidelines avail- 
able concerning the spectral characteristics of a stellar-mass black 
hole accreting at such super-Eddington rates. They could quite 
conceivably mimic the spectrum of a more massive bla<;k hole 
accreting at sub- Eddington rates. 
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2 BINARY EVOLUTION CALCULATIONS 

In our previous study related to ULXs (Podsiadlowski, Rap- 
paport, & Han 2003; hereafter PRH) we developed two sep- 
arate components of the calculations necessary for a binary 
population synthesis of these objects. In the first, we started 
with a very large set of massive primordial binaries and gen- 
erated a much smaller subset of these that evolved to contain 
a black hole and relatively unevolved companion star. The 
product was a set of "incipient" black-hole X-ray binaries 
with a particular distribution of orbital periods, Porb, donor 
masses, M2, and black-hole masses, Mbh, for each of a num- 
ber of different sets of input assumptions (sec, e.g.. Fig. 2 of 
PRH). For this part of the calculation, we employed various 
"prescriptions" , based on single star evolution models for the 
primary, simple orbital dynamics associated with wind mass 
loss and transfer, assumptions about the magnitude of the 
wind mass loss from the primary as well as from the core of 
the primary after the common envelope, and natal kicks dur- 
ing the core collapse and formation of the black hole. Simple 
energetic arguments were used to yield the final-to- initial or- 
bital separation during the common-envelope phase wherein 
the envelope of the primary is ejected. Here we utilized a pa- 
rameter A, which is the inverse of the binding energy of the 
primary envelope at the onset of the common-envelope phase 
in units of GMiMe/ Ri, where Mi, and R\ arc the total 
mass, envelope mass, and radius of the primary, respectively. 
This parameter strongly affects the final orbital separation 
after the common-envelope phase, where smaller values of 
A correspond to more tightly bound envelopes, and hence 
more compact post-common-envelope orbits. 

Conventional energetic arguments for the ejection of the 
common envelope yield the following expression for initial- 
final orbital separation: 

(^) -^(^^ + ^)" ' (2) 
Vai/cB Ml \ r7cBArL/ 

(e.g., Webbink 1985; Dewi & Tauris 2000; Pfahl, Rappa- 
port, & Podsiadlowksi 2003), where the subscripts "1", "c", 
and "e" stand for the progenitor of the black hole, its core, 
and its envelope, respectively, and "2" is for the progenitor 
of the "donor star" in the black-hole system. The quantity 
ri, is the Roche lobe radius of the black-hole progenitor in 
units of Oj, ?7cE is the fraction of the gravitational binding 
energy between the secondary and the core of the black-hole 
progenitor that is used to eject the common envelope, and 
A is defined above. For typically adopted parameter values, 
A ~ 0.01 - 1 (e.g., Dewi & Tauris 2000), t^ce ^ 1, and 
tl — 0.45 — 0.6 (for an assumed mass ratio between the 
black-hole progenitor and the companion in the range of 
■~ 2 : 1 ^ 15 : 1), the second term within the parentheses 
in eq. (2) dominates over the first. In this case, we find the 
following simplified expression for af/ai: 

f^) -^fT^)^^A=. 0.005 f^')A, (3) 

where the leading factor is r_L/2 ~ 1/4, while the factor in 
parentheses involving the black-hole progenitor is ~ 0.020 ± 
0.002 Mg^ for virtually all of the progenitors we consider. 
This explains why the large majority of the incipient black- 



hole binaries found by PRH resulted from an initially very 
wide orbit (Porb ~ years — when the primary attains radii 
of ~ 1000—2300 Rq) preceding the common-envelope phase 
in order to avoid a merger between the secondary and the 
core of the primary. 

The final orbital separations (and corresponding orbital 
periods) after the common envelope depend on the initial 
orbital separations according to eq. (3). These, in turn, de- 
pend upon a number of issues associated with the evolution 
of massive stars and their wind loss characteristics which 
were explored in PRH. We expect, however, that the dis- 
tribution of incipient black-hole binaries in the Porb — M2 
plane will lie largely within an envelope bounded at the top 
by Porb.t oc {aiA42^)^^^ (motivated by eq. [3]), and at the 
bottom by values of Porb.b that just barely avoid a merger 
between the core of the black-hole progenitor and the com- 
panion star. From the results of PRH (see, e.g., their Figs. 
2 &: 3) we can fit the following semi-empirical expressions 
to the regions in the Porb — M2 plane where the incipient 
black-hole binaries are located: 

for M2 > 1 Mq. The expression for the lower limit on Porb 
results from the radius-mass relation for stars on the main 
sequence and the functional dependence of the Roche-lobe 
radius on mass ratio. Furthermore, we find empirically that 
the lower limit on the donor masses in these incipient sys- 
tems is given by: 

M2,min ~ ^ Mq . (6) 

This expression results from the fact that the product of 
M2A must exceed a certain value (see eq. [3]) in order to 
avoid a merger. A more extensive discussion of the genera- 
tion of the incipient black-hole X-ray binaries and the rel- 
evant references are given in PRH. For other references on 
this subject see, e.g.. King & Kolb (1999); Fryer & Kalogera 
(2001); Portegies Zwart et al. (1997); and Nelemans & van 
den Heuvel (2001). 

For the second part of the PRH study we calculated a 
small but unique grid of 19 black-hole binary models evolv- 
ing through the X-ray phase, i.e., where the donor star trans- 
fers mass to the black hole. These calculations are done with 
a full Henyey stellar evolution code so that at least the be- 
havior of the donor star throughout its mass-transfer phase 
should be quite accurate and realistic. Nonetheless, there are 
still a number of uncertainties which involve such issues as 
whether the Eddington limit is strictly adhered to, and when 
the mass transfer onto the black hole becomes "transient" 
in nature; we investigate a number of these issues in this 
work. In PRH we systematically explored a subset of case 
A binary evolution models where the donor star is initially 
unevolved at the time when mass transfer onto the black 
hole commences. Fourteen of the PRH models were for this 
early case A binary evolution. A few additional models were 
computed for cases where hydrogen had been partially de- 
pleted (also case A) or completely exhausted (early case B) 
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Figure 1. Potential X-ray luminosities as a function of time for 
14 black-hole X-ray binary evolution sequences; each color cor- 
responds to a different initial donor star mass ranging from 2 to 
17 Mq (adapted from PRH). The trend in the evolutionary se- 
quences from higher to lower overall X-ray luminosity corresponds 
to decreasing initial donor mass. In all cases the donor star is un- 
evolved at the start of mass transfer, and the black-hole mass is 
IOMq. The spiky feature at the end of each evolution corresponds 
to the donor star ascending the giant branch. The duration of this 
latter phase is ~5% of the entire evolution. 
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Figure 2. Initial donor masses and orbital periods for the 52 
black-hole binary evolution models. 



when mass transfer commenced. PRH found that the X-ray 
lifetime in these latter systems is 1 to 2 orders of magnitudes 
shorter than for systems experiencing case A mass transfer 
(see Table 1 of PRH and Figs.l3land lll|l . Moreover, for large 
initial mass ratios, mass transfer would be dynamically un- 
stable. We therefore expect that these more evolved systems 
(cases B and C) will make only a relatively small contribu- 
tion to the overall population of ULXs. 



We reproduce in Fig. plots of the "potential" X-ray 
luminosity (L§°*) for 14 binary evolution sequences with ini- 
tial donor masses of 2 through 17 Mq from the original 
models presented in PRH (for the case of initially unevolved 
donors). The "potential" X-ray luminosity refers to the en- 
ergy output expected in the absence of the Eddington limit; 
i.e., if the X-ray luminosity^, L^, were limited only by the 
mass transfer rate, M, and the energy conversion efficiency 
which is dictated by the instantaneous spin of the black hole 
(Bardeen 1970). We assumed that each black hole starts with 
a mass equal to 10 Mq and with zero spin (i.e., j = 0, where 
< J < 1 is the dimensionless spin parameter). We reiter- 
ate that in each model shown in Fig. the donor star is 
assumed to start on the main sequence when mass transfer 
commences. We also note that these binary models include, 
in addition to M and Lx, the evolution of Porb, A/2, Mbh, 
and the spin of the accreting black hole. 

In the present work we have significantly augmented our 
grid of black hole binary evolution models from 19 to 52, and 
in the process consider case A mass transfer in a more sys- 
tematic way. The additional models follow the same mass 
grid as the original 19, except that for each mass we now 
allow for three additional evolutionary states of the donor 
star at the onset of mass transfer. This effectively takes into 
account a significant portion of the expected orbital period 
distribution among the incipient systems. The different ini- 
tial evolutionary states of the donor are characterized by the 
hydrogen mass fraction in the core, Xc- Thus, for each mass, 
the models are for values of Xc ~ 0.7, 0.35, 0.2, and 0.1. The 
corresponding orbital periods at the start of mass transfer 
range from about 0.6 — 4 days. The location of the starting 
points for these models in the M2 — Porb plane are shown in 
Fig. 121 To our knowledge, these are the only properly com- 
puted evolution models for black-hole binaries during their 
X-ray phase. 

Plots of vs. time since the birth of the incipi- 

ent binary for a sample of our new evolutionary sequences 
are shown in Fig. |3 these are all for the case where the 
initial donor mass is 10 Mq, but the initial H fraction 
is different; the black, red, blue, and green curves are for 
Xc — 0.7,0.35,0.2, and 0.1, respectively. Note that, as 
the donor star is progressively hydrogen depleted at the 
start of mass transfer, the duration of the mass transfer 
phase becomes systematically shorter, while the peak values 
of become correspondingly higher. The systematically 
later times of the onset of mass transfer simply reflect the 
time that the donor star requires to evolve and fill its Roche 
lobe. The duration of the main mass-transfer phase becomes 
shorter, both because the nuclear burning timescale which 
sets the rate of evolution of the star becomes shorter (and 
hence leads to a higher mass-transfer rate) and because the 
remaining core hydrogen-burning lifetime decreases, since 
less H is available. The rest of the 52 evolutionary sequences 
are not shown here in the interest of space, but are used in 



^ Here, as throughout the paper, all calculated luminosities refer 
to the total accretion luminosity without regard to the region 
of the X-ray band in which it emerges. Thus, the 2 — 10 keV 
luminosity would likely be somewhat lower than the values cited; 
see §4 for a discussion of this issue. 
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Figure 3. Potential X-ray luminosities as a function of time for 
5 black-hole X-ray binary evolution sequences; each color corre- 
sponds to a different evolutionary state of the donor star when 
mass transfer commences. In all cases the initial masses of the 
donor and black hole are 10 Mq and 10 Mq , respectively. The 
black, red, green, and blue curves correspond to the central hy- 
drogen abundance of the donor star at the onset of mass transfer 
of 0.7, 0.35, 0.2, and 0.1, respectively. For comparison, the dot- 
ted curve shows the potential X-ray luminosity for a 10 Mq sec- 
ondary that has already evolved off the main sequence (so-called 
early case B mass transfer), where mass transfer is driven by the 
expansion of the secondary as it crosses the Hertzsprung gap on 
a thermal timescale. The duration of the mass-transfer phase is 
not resolved in the figure, as it lasts only ~ 50, 000 yr. 



the population synthesis study presented in this work, and 
are available in digital form upon request to the authors. 

Experience has shown that population synthesis calcu- 
lations for systems containing collapsed stars typically re- 
quire > 10'' systems in order to achieve results of high statis- 
tical quality (see, e.g., Howell, Nelson, & Rappaport 2002). 
Since each of the 52 models we have computed required some 
half hour of cpu time — as well as some hand holding — it is 
therefore impractical at this time to run thousands of black- 
hole binary models with a Henyey code to describe the pop- 
ulation. However, we have found from our models that the 
evolutions in form a nearly self-similar set. That is, if we 
take any two of the evolutionary models, but especially ones 
that are close in initial AI2 and Porb, then the plots of 
vs. evolution time can be scaled in time and in so they 
are nearly the same (see, e.g., Figs. and 1^. This is quite 
different from the case of binary evolution for systems con- 
taining neutron stars (see, e.g., Podsiadlowski, Rappaport, 
& Pfahl 2002) where the evolution tracks for initially similar 
systems can diverge dramatically. Therefore, we have made 
use of this self-similar behavior to develop an interpolation 
scheme to produce an effectively much larger set of models. 
Thus, if we wish to know the evolution of a binary system 
starting with values of Porb and M2 that are within the grid 
boundaries, but are not located at one of the 52 "nodes" , we 



carry out the following variant of a bi-linear interpolation. 
First, a weighting factor for each of the four nearest grid 
models (in the Porb — M2 plane) is computed by linearly in- 
terpolating in both Porb and M2. Once this has been done, 
the logarithms of both the beginning and the end times {ti 
and tf) of the four evolutions are interpolated using the 
weighting factors; this then provides values of < ti > and 

< tf > to use in the interpolated evolution. Each of the four 
nearest evolutionary models then has its time axis stretched 
and shifted so that it begins and ends at < ti > and < tf >, 
respectively. For each time between < ti > and < tf > we 
then do a logarithmic weighting of the quantity of interest 
(e.g., Lx) for the four nearest models, using the same weight- 
ing factors as described above. We have visually inspected a 
substantial number of the interpolated evolutions as a check 
that they behave sensibly. 

The black-hole binaries that we produce may account 
for a significant fraction of the most luminous X-ray sources 
born in active star forming regions. However, there are sev- 
eral classes of objects not generated in our population syn- 
thesis. Not included are the compact black-hole binaries, 
with low-mass companions, that are transient sources and 
are observed in abundance in our own Galaxy (see the dis- 
cussion in PRH). We argue below that these, in fact, may be 
a separate population and not what is being observed in star 
formation regions. Moreover, observations of these sources 
in our own Galaxy indicate that they tend to be limited to 

< 10''* ergs , so they may not contribute substantially to 
the ULX population. We do not consider the evolutionary 
phase before Roche-lobe overflow when the black hole is fed 
by accretion from the stellar wind of the companion. Such 
systems are likely to be at the low end of the luminosity 
functions we generate in this work, e.g., ~ 1O'''*(Mbh/M0)'^ 
ergs s~^, where Mbh is the mass of the wind-accreting black 
hole. Also, we do not produce systems where the black hole 
is formed in very wide primordial binaries without a com- 
mon envelope being involved, as in the Voss & Tauris (2003) 
scenario. If objects collapsing to form black holes have natal 
kicks, and if these kicks are appropriately scaled down by 
the mass of the collapsed star from the case of neutron star 
kicks (e.g., to preserve a constant recoil momentum), then 
the subsequent orbit could fortuitously remain bound and 
become highly eccentric. Eventually, the donor star would 
evolve to fill its Roche lobe and mass transfer would then 
commence. Such systems will, however, tend to have much 
longer orbital periods than the ones we generate. This has 
two consequences: (1) many of these will result in unstable 
mass transfer if the donor star has a mass exceeding that 
of the black hole; and (2) the mass-transfer phase will be 
extremely short lived. Finally, we are not considering any 
systems containing neutron stars since these are not likely 
to be able to explain the most luminous ULXs. On the other 
hand, such systems will surely contribute to the luminosity 
function up to at least ~ 5 x 10^* ergs s"'^. 
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Figure 4. Simulated evolution of the X-ray luminosity function of black hole binaries with time since an impulsive star formation event. 
"Luminosity" refers to the total accretion luminosity without regard to the region of the X-ray band in which it emerges. Thus, the 
2 — 10 keV luminosity would likely be somewhat lower than the values indicated in the plots. Evolution tracks from 10^ X-ray sources 
were computed and then registered in each of the 700 X 700 pixels that are traversed. The colors crudely represent the logarithm of 
the relative populations, with purple through red corresponding to ratios of 200 to 1. The mass of the black hole at the start of mass 
transfer is taken to be 10 Mq. Panel a: Model E. The appropriate Eddington limit of each system has been applied. The upper red 
bar-like feature between logt ~ 7.2 — 8.2 represents systems with donor stars on the giant branch with higher mass transfer rates; the 
larger black- hole masses later in the evolution and higher He fraction raise i^Edd- Panel b: Model B. All specifications are the same as 
for panel (a), except that the Eddington limit is allowed to be violated by up to a factor of 10 (see text for details), and each system 
radiates the luminosity that corresponds to the mass transfer rate and the energy conversion efficiency according to the instantaneous 
spin of the black hole (Bardeen 1970). As in panel (a), the red features at the higher Lx values (at logtev^7.2) are systems with the 
donor on the giant branch. Panel c: Model S. All specifications are the same as for panel (b), except that the initial mass of the donor 
star is chosen from a Salpeter initial mass function (Salpeter 1955). This approximates the output from our binary population synthesis 
code for the case where A ~ 0.5. Note that compared to panels (a) & (b) there are many fewer high Lx systems due to the emphasis on 
the lower-mass donors. Panel d: Model TS. All specifications are the same as for panel (c), except that transient behavior of the X-ray 
sources is approximately taken into account (see text for details). Such systems are assumed to be in an "on" state only 3% of the time, 
but with a value of M that is 30 times higher than would be the case if the X-ray emission were steady. Note that most of the systems 
with Lx < 10'^^ ergs s~^ no longer appear on the figure, but have been systematically shifted to higher Lx, producing the yellow ridge 
for main-sequence donors. The faint red feature at times exceeding lO'' yr with high Lx, corresponds to transients with the donor stars 
on the giant branch. 
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3 POPULATION SYNTHESIS 
CALCULATIONS 

We now use the tools described above to generate large 
populations of black-hole binaries of the type likely to be 
produced in active star formation regions. The first step is 
to decide which of the outputs to use from the code that 
generates the incipient black-hole binaries. Since there are 
many uncertainties that go into these calculations, we have 
decided to use somewhat more general distributions of in- 
cipient systems in the Porb — M2 plane that are inspired by 
the results we found in PRH, rather than taken directly from 
the specific output for any particular model. We utilize three 
generic types of distributions for the incipient black-hole bi- 
naries. In the first, we take Porb to be uniformly distributed 
between the two curves given by eqs. (1) and (2), and M2 to 
be uniformly distributed between 2 and 17 Mq. This allows 
for the broadest contributions from all the binary evolution 
models that we have calculated, and would directly imply 
a large value of the parameter A. The second case we con- 
sider is one where Porb is distributed uniformly as above, 
but M2 is distributed in the same way as the stellar initial 
mass function as deduced by, e.g., Salpeter (1955), Miller & 
Scalo (1979), and Kroupa, Tout, & Gilmore (1993). This is 
a very good approximation to what we find for the high A 
case, e.g., A > 0.5 (PRH). In this case the envelope of the pri- 
mary is relatively easy to eject and secondaries over a wide 
range in mass are able to successfully get past this phase. 
Their mass distribution then also roughly resembles that of 
the primaries since we utilize a flat mass ratio distribution 
(sec PRH). The third case we consider is one where A has 
a small value, e.g., 0.08. In this case only the more mas- 
sive secondary stars are successful in ejecting the envelope 
of the primary. For the specific value of A = 0.08 the incip- 
ient population of black-hole binaries has secondary masses 
largely confined to A42 }^ 6 A4q , and the masses arc roughly 
uniformly distributed above this value (see Fig. 2 of PRH) . 

Throughout this work we have somewhat arbitrarily 
adopted a fixed black-hole binary production rate of TZbh — 
10~^ yr~^ and cite a representative value of A ~ 0.1 in or- 
der to normalize all of our luminosity functions. In fact, our 
earlier binary population synthesis study (hereafter "BPS" ; 
PRH) yielded the following values: TZbh — 10~^ yr~^ and 
A ~ 0.08 for an incipient black-hole binary population 
characterized approximately by donor masses uniformly dis- 
tributed between 6 and 17 Mq, and 7?.bh — 3 x 10^® yr^^ 
and A ~ 0.5 for the case where the donor masses are approx- 
imately distributed according to the Salpeter (1955) distri- 
bution over the mass range 2 — 17 Mq. On the other hand, 
no particular BPS model output corresponded closely to our 
current choice of a uniform distribution of donor masses be- 
tween 2 and 17 Mq, and thus we cannot cite BPS values 
for 7?^BH and A for this model. In all of our BPS models 
(PRH) the core-collapse supernova rate was taken to be 
TtsN = 10^^ yr^^. We comment later in the text on how 
appropriate we think our adopted rate of 10~^ yr~^ is for 

7?-BH- 

In addition to the three difi^erent distributions of incip- 
ient black hole binaries in the Porb — M2 plane described 
above, there are two different assumptions we have made 



regarding the maximum X-ray luminosities that can be ra- 
diated by a black-hole binary. In the first of these, the max- 
imum value of Lx is taken to be just that given by the Ed- 
dington limit which, in turn, is governed by the mass of the 
black hole and the hydrogen/helium composition of the ac- 
creted material. In the second, we allow for the possibility 
that luminosities up to 10 times the nominal value of LEdd 
can be attained (Begelman 2002; Ruszkowski & Begelman 
2003). Finally, we consider one additional case where the lu- 
minosity of the X-ray sources is, under certain physical con- 
ditions, transient in its behavior. The transient behavior is 
thought to arise from the well-known thermal-ionization disk 
instability (Cannizzo, Ghosh, & Wheeler 1982; van Paradijs 
1996; King, Kolb, & Burderi 1996; Dubus, Hameury, & La- 
sota 2001; Lasota 2001). In particular, if the X radiation 
from the central source is not able to maintain a sufficiently 
high temperature for the accretion disk (especially for the 
outer portions), the mass transport through the disk would 
be unstable, and vice versa. In general, this affects mostly 
the lower luminosity sources with L^, < 10^^ ergs s~^. The 
specific prescription we used to determine whether a partic- 
ular model at a given time would exhibit transient behavior 
is described in detail in PRH (sec also §4). When one of 
our sources is in an evolutionary phase where it would be 
a transient we simply allow, in a somewhat ad hoc maimer, 
the mass transfer rate onto the black hole to increase by 
a factor of 30 above the long-term average rate of matter 
flowing into the disk, but then give that source a probabil- 
ity weighting of only 1/30 in the population, indicating that 
it would be in an "on" state only ~3% of the time. More 
empirically correct factors to represent transients might be 
as low as 0.01 for the duty cycle (McClintock &: Remillard 
2004; and references therein) and up to 100 times the long- 
term average mass transfer rate. However, as we shall show, 
transient systems play a relatively minor role in our models 
of the ULX population, and that would not change even if 
we had adopted these latter parameters. 

In summary, the six models that we have chosen to work 
with are deflned as follows: 

• Model E is defined to have a uniform incipient population 
in the Porb — M2 plane, and the luminosity during the X-ray 

phase is limited to LecIci; 

• Model B is the same, except with a "Begelman" luminos- 
ity enhancement such that Lx is permitted to be as large as 
10 Lfidd; 

• Model TB is the same as Model B, but transient source 

behavior is also included; 

• Model S is the same as Model B except that M2 is dis- 
tributed according to a Salpeter (1955) mass function rather 
than uniformly in mass; 

• Model TS is the same as Model S, but transient source 

behavior is also included; 

• Model 6 is the same as Model B, except that only values 
of M2 > 6 are included. 

For each of these models we generate a luminosity func- 
tion vs. time, tcv, since an impulsive star formation event. 
We choose, via Monte Carlo draws, 10'' incipient black-hole 
binaries in the Porb — M2 plane. To the extent that a par- 
ticular value of Porb and M2 lies within our grid of models, 
we compute an interpolated evolution from the 4 nearest 
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Table 1. Summary Model Parameters 



Name 


Model 


N{> lO^^ergss-i) 


N{> lO^^-S ergss-i) 


r{«) 
-^x.tot 


-^x,tot 


E 


Eddington limited luminosity 


30 


0.2 


4 X 10*0 


1.5 X 10*1 


B 


Begelman super Eddington factor 


30 


6 


6 X 10*° 


2.3 X 10*1 


TB 


Transient behavior included 


30 


6 


6 X 10"° 


2.3 X 10*1 


S 


Salpeter IMF 


4 


0.7 


1.0 X 10*° 


1.5 X 10*1 


TS 


Transient behavior included 


4 


0.7 


1.0 X 10*0 


1.5 X 10*1 


6 


initial donor masses > 6Mq 


40 


8 


9 X 10*0 


2.6 X 10*1 



(a) Lx.tot is in units of ergs s~i and represents the maximum in the total X-ray luminosity of a galaxy (due to BH-binaries) after a 
starburst event yielding 10^ core-collapse SNe. (b) Lx.tot is in units of ergs s~l and represents the steady-state total X-ray luminosity 

of a galaxy (due to BH-binaries) undergoing 0.1 core-collapse SNe per year. 
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Figure 5. Calculated X-ray luminosity function of black- hole bi- 
naries at five different epochs after an impulsive star formation 
event (units are BH-binaries per luminosity bin). The time labels 
on each curve indicate log(tcv) in units of years. Solid curves are 
for Model B (see Fig. IH-panel b). The two dashed curves are for 
Model E (i.e., Eddington limited X-ray luminosities; see Fig. IH- 
panel a) at log(tev) = 7.0 and 7.5; the curves for all later times 
are nearly indistinguishable from those of Model B. The distribu- 
tions are normalized to a star formation event which yielded 10® 
core-collapse SNe. 
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Figure 6. Calculated cumulative X-ray luminosity function of 
black-hole binaries at five different epochs after an impulsive star 
formation event. The time label on each curve indicates log(tev) 
in units of years. Solid curves are for Model B (see Fig.Bl-panel a). 
The two dashed curves are for Model E (i.e., Eddington limited 
X-ray luminosities; see Fig. IH-panel a) at log(icv) = 7.0 and 7.5; 
the curves for all later times are nearly indistinguishable from 
those of Model B. The normalization of the curves is the same as 
in Fig. 121 



models, as described in detail in §2. To store the results, 
we set up an array in the log — log iev plane, containing 
700 X 700 elements covering the range 10^'* < < 10*1 ergs 
s"^ and 10'' < tcv < 10^'^ yr in equal logarithmic bins. Each 
time that one of our evolution tracks crosses an element of 
this array, a value of unity is added to that element. After 
this operation has been completed for all 10^ interpolated 
evolution tracks, the result is displayed as an image repre- 
senting the unnormalized evolving luminosity function for 
each model. These are shown for Models E, B, S, and TS, 
in Fig. |1| The luminosity function for Models 6 and TB are 
not shown in the interest of space. 

Model E (Fig.|l] panel (a)) shows a large concentration 
of highly luminous sources (^10'^^ ergs s~i; i.e., ULXs) at 
10 — 30 Myr. These result from systems with the initially 
more massive companion stars. Between ~30 and 300 Myr, 
the peak in the luminosity function falls monotonically from 
~ 10^^ to < lO^'^ ergs s~i as systems with the more massive 



donors run their evolutionary course, leaving longer-lived, 
lower-luminosity systems with initially lower-mass donors. 
Numerous systems with L-x exceeding Z/Edd for a neutron 
star are still present after 100 Myr. Note the red group of 
systems at luminosities above the continuous distribution 
starting at ~ 15 Myr and continuing to > 1 Gyr; these result 
from the relatively short-lived phase when the donor stars 
(of all initial masses) ascend the giant branch and produce 
mass transfer rates that increase by 1 - 2 orders of magnitude 
over the preceding phase of the binary evolution. Finally, 
note that the apparent upper cutoff in the luminosity func- 
tion seems to increase by about a factor of 2 for evolution 
times later than ~ 20 Myr — this results from the general in- 
crease in the black-hole mass as they grow by accretion, and 
the changing chemical composition of the accreted material 
which is becoming He enriched. 

The evolving luminosity function for Model B (Fig. 01 
panel (b)) is very similar to that of Model E, except for 
the pronounced increase in the allowed upper limit to Lx, 
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Figure 7. Number of luminous black-hole X-ray binaries as a 
function of time after an impulsive star formation event. Solid 
curves represent the number with > 10^® ergs s~^; dashed 
curves are for Lx > 3 X 10^^ ergs s"'^. The normalization of the 
curves is the same as in Fig. |21 The peak value for each curve is 
listed in Table 1. 




Log (Time - yr) 

Figure 8. Total luminosity from black-hole X-ray binaries as 
a function of time after an impulsive star formation event. Re- 
sults for Models E, B, S, TS, and 6 are shown. The distributions 
are normalized to a star formation event that yielded 10^ core- 
collapse SNe. The peak value for each curve is listed in Table 1. 



by a factor of 10 - as defined by the model. This model 
allows for a significant population of sources with Lx up to 
~ 3 X 10'^^ ergs s~^ and a tailing off population to values of 
Lx — 2 X 10*° ergs s~^ for times up to 30 Myr after a star 
formation event. We note that the luminosity function for 
Model B (Fig. 2] panel (b)) extends well past the peak in 
the distribution at all times. Thus, increasing the "Begelman 
factor" by more than a factor of 10 would not substantially 
enhance the population of ULXs. 

In Fig. 2] panel (c) we present the evolving luminosity 
function for Model S, where the initial donor star masses are 
distributed according to a Salpeter IMF, i.e., with a steeply 
decreasing population with increasing mass. In such a sce- 
nario, the X-ray luminosity function is much less strongly 
peaked at times ^30 Myr, i.e., not dominated by a nar- 
row range of luminosities, than in the previous models. In 
particular, there are relatively fewer candidate ULX sources. 
This is a natural consequence of the smaller numbers of mas- 
sive donor stars. Finally, the evolving luminosity function for 
Model TS is presented in Fig. 0] panel (d). In this model, 
the Salpeter mass distribution is again applied to the in- 
cipient donor stars, but transient source behavior due to 
the thermal/viscous disk instability is now included (in an 
approximate way) . Note that most of the lower luminosity 
sources (i.e., with Lx ^ 10^*^ ergs s~^) are missing from the 
evolving luminosity function. These have been boosted up 
by a factor of 30 in Lx, but are given a weighting in the 
"image" of only 1/30. The effect of the transients is most 
pronounced for times > 200 Myr where they appear as a 
yellow ridge in the Lx — tcv plane with Lx in the range of 
1q38_]^q39 gj.gg g-i f,j-Q^g g^igQ ^jjg faint red feature appearing 

at times > 10® yr which extends into the ULX luminosity 
range. These are due to donor stars of initially lower mass 
that are on the giant branch and exhibit transient behavior 
(GRS 1915-1-105 in our Galaxy may be an example of such 
a system). 



The evolving luminosity functions presented as images 
are visually informative and interesting, but a somewhat 
more quantitative way to view these distributions is as one- 
dimensional plots of either the luminosity functions at spe- 
cific times or as the evolution with time of the number of 
sources exceeding a given Lx. The first of these is shown in 
Fig-El Here, the luminosity function is plotted for Model B 
at 6 different epochs after a star formation event; these are 
spaced logarithmically every factor of vTO in time, starting 
at 10^ yr. Also shown for comparison as a set of dashed lines 
are luminosity functions for Model E at log(fcv) ~ 7.0 and 
7.5; the curves for all later times are nearly indistinguishable 
from those of Model B. The curves are normalized to 10^ 
core-collapse supernovae associated with the star formation 
event (including those in binary as well as single stars), and 
to a ratio of the black-hole production rate to that for core- 
collapse SNe, 7?.bh/7?.sn, equal to 10"*. The most prominent 
appearance of the ULXs, as could be seen qualitatively in the 
luminosity function images, occurs between 10 and 30 Myr. 
The same luminosity functions, but plotted as cumulative 
distributions, are shown in Fig. |^ From this plot it is easy 
to see that between 10 and 30 Myr a galaxy, undergoing suf- 
ficient star formation to produce 10® SNe impulsively (e.g., 
within < 10 Myr) , could be expected to harbor more than 
a dozen ULXs containing stellar-mass black holes. Also, we 
show in Fig. □ the number of ULXs ( > lO^'' and > 10^®'^ 
ergs s~^) vs. time since a star formation event for Models 
B, E, S, TS, & 6. These curves have also been normalized to 
10^ SNe and to 7?.bh/7?.sn = 10~*. Finally in this regard, we 
show in Fig. |H|the integrated X-ray luminosity correspond- 
ing to the source numbers given in Fig. Q Except for Model 
S, typical integrated X-ray luminosities of 4 — 9 x 10*" ergs 
s~^ are attained following a large star formation episode 
(normalized to lO'^ core-collapse SNe). 

The type of information contained in the 'data' files 
used to produce Fig. 21 can also be used to construct the 
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Log (Luminosity - ergs/sec) 

Figure 9. Luminosity functions for a galaxy undergoing contin- 
uous star formation at a fixed age ^ 10* yr. The core-collapse 
SN rate is taken to be 0.01 yr~^. The results are derived from 
Fig. |1] by weighting each column of pixels by the time interval 
represented by that pixel, and then summing each row of pixels 
over all columns. Models B, E, S, TS, and 6 are defined in the 
text. 



black-hole binary X-ray luminosity function for any history 
of star formation. To illustrate this, we compute the lumi- 
nosity function for a galaxy undergoing continuous, uniform 
star formation for times exceeding billions of years. To do 
this we merely sum over the columns of the — tcv ma- 
trix, multiplying each column by its width in time (noting 
that the widths of the columns increase logarithmically with 
tev). The results are shown in Fig. |3|for an assumed forma- 
tion rate of black-hole binaries of 7?.bh = 10~^ yr~^ (PI^H), 
which in turn utilizes a SN rate of 0.01 yr~^ (e.g., Cappel- 
laro, Evans, & Turatto 1999). The results in Fig.El for three 
of our models (B,E, and 6), indicate that a typical normal 
spiral galaxy such as our own would have about a half dozen 
ULXs at any given time. Models S and TS, with a smaller 
population of initially massive stars, produce of order one 
ULX, in closer agreement with the observed number in our 
Galaxy. 

We can also use the results of Fig.|5]to compute the in- 
tegrated X-ray luminosity in the continuous star formation 
scenario for each of our models. We list the results in Ta- 
ble 1 — scaled up by a factor of 10 for better comparison to 
typical active star forming galaxies with an assumed core- 
collapse SN rate of 0.1 yr~^. It is instructive to compare our 
results with the empirical findings of Gilfanov, Grimm, & 
Sunyaev (2004) and RanaUi, Comastri, & Seti (2003) who 
cite their findings in terms of the relation between the mas- 
sive star formation rate (SFR), expressed in units of Mq 
yr~^ for stars with M > 5 Mq, and the total X-ray luminos- 
ity Z/x,tot. To intercompare our theoretically derived results 
with their empirical findings, it is helpful to have a conver- 
sion factor between the rate of core-collapse SNe (7?.sn) and 
the SFR, which we estimate to be 7^sN ^ 0.025 SFR for 
a Salpeter-type IMF. With this conversion we find that our 
adopted value of TZsn ~ 0.1 yr~^ corresponds to an SFR of 
~ 4 Mq yr~^. For this particular value of SFR Gilfanov et 
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Figure 10. Simulated representation of the luminous black-hole 
X-ray binaries in a galaxy where a star formation wave was trig- 
gered by a catastrophic event at the center some 5 X 10* years 
ago. The top and bottom panels were produced for Models B 
and 6, respectively. The wave of star formation is assumed to 
propagate outward at a constant speed. The surface density of 
interstellar gas used to form the stars is taken to be a constant 
per unit surface area. The X-ray luminosities were chosen via a 
Monte Carlo technique from the distribution for Model B (Fig. 
l4l-Danel b) and Model 6. This spatial distribution of sources is 
reminiscent of Chandra images of the Cartwheel galaxy. 



al. (2004) and RanaUi et al. (2003) find an average value of 
ix.tot — 2 X 10'"' ergs s^^. Our tabulated values of Lx,tot 
(Table 1) for continuous star formation are typically an or- 
der of magnitude larger than than the empirical value. The 
values for Lx,tot in Table 1 for impulsive star formation are 
in better agreement with those of Gilfanov et al. (2004) and 
RanaUi et al. (2003), but there our normalization of 10® 
core-collapse SNe, which corresponds to, e.g., 7?.sn = 0.1 
yr~^ for a finite time of 10^ years, probably underestimates 
the actual number of SNe contributing to the observed pop- 
ulation of luminous X-ray sources. On the other hand, the 
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assumption of a completely continuous star formation sce- 
nario likely overestimates the number of luminous systems. 
On balance, we conclude that the comparison of our results 
(Table 1) to those of Gilfanov et al. (2004) and Ranalli et al. 
(2003) indicates that our adopted formation rate of black- 
hole binaries of 10~® yr~^ for 7?.sn = 0.01 yr~^ is likely too 
high by a factor of ~ 3. This probably also indicates that our 
normalization to an effective value of A = 0.1 (the binding- 
energy parameter of the BH progenitor) is likely a bit too 
large; perhaps a better selection would be A = 0.06 — 0.07. 

Finally, we illustrate how, for external galaxies, our re- 
sults can be used to relate the spatial offsets between ULXs 
and current star forming regions. For this we use the evolving 
luminosity functions after an impulsive star formation event 
(e.g.. Fig. 0J. In this exercise we make a crude model for 
the annular concentration of luminous X-ray sources in the 
Cartwheel galaxy. We envision an unperturbed spiral galaxy 
with a uniform surface column density of interstellar gas as 
the initial condition. We assume that some 5 x 10* years ago 
a smaller galaxy passed through the center of the spiral and 
triggered a wave of star formation (Hernquist & Weil 1993) 
which has just now reached a radius marked by numerous 
HII regions and luminous X-ray sources. We choose random 
locations over the surface of the disk (per unit area); the 
radial distance of each location then specifies the time since 
the star formation wave has passed. The results shown in 
Fig. m panel (b) are then used as a probability distribution 
to choose an X-ray luminosity. This is repeated until 100 X- 
ray sources have been chosen for the realization. The results 
are shown in the top panel of Fig. 1101 The yellow ring marks 
the outer 10% in radius. Note that while the luminous X-ray 
sources are indeed concentrated in the outer ring, there are 
still a number of luminous sources (>10^** ergs s"^) well 
inside the ring. We then repeated this exercise for Model 
6 where initially more massive donor stars are emphasized; 
for this model the sources are even more concentrated in the 
outer 10% of the simulated galaxy (lower panel in Fig. llO^ . 



4 DISCUSSION AND CONCLUSIONS 

We have computed a unique grid of 52 black-hole binary evo- 
lution models covering donors in the mass range of 2 — 17 Mq 
and 4 different evolutionary phases at the onset of mass 
transfer from the ZAMS to the TAMS. We chose the incipi- 
ent black-hole binaries from regions of the Porb — M2 plane 
that were suggestive of the results of the binary population 
synthesis study by PRH. Several different assumptions were 
made in choosing the donor masses in the incipient popula- 
tion of black-hole binaries (i.e., Model sets {E,B,TB}, {6}, 
and {S,TS}). For each incipient binary chosen, we utilized an 
evolution track that was interpolated from among the binary 
evolution models in our grid. In this way we computed lumi- 
nosity functions for large sets of simulated black-hole bina- 
ries. The masses of all the black holes in our evolutionary cal- 
culations were taken to be initially 10 Mq. In future studies 
the black-hole mass should be taken directly from the pop- 
ulation synthesis calculations, and would then range from 
~ 5 — 15 Mq (based on observational constraints; e.g., Mc- 
Clintock & Remillard 2004). The lower black-hole masses. 



when paired with the higher-mass donors (e.g., ^12 Mq), 
would lead to dynamically unstable mass transfer and would 
then not enter the population. The inclusion of lower-mass 
black holes will have some modest effect, but certainly not a 
dominant one, on the luminosity functions that we calculate 
(see also Fig. II II and the associated discussion). 

Our luminosity functions are sensitive to two princi- 
pal input parameters. First, the value of A (related to the 
binding energy of the progenitor of the black hole) directly 
affects the mass distribution of donor stars in the black-hole 
binaries. As discussed in the text, there are significant un- 
certainties in the appropriate values of A to use. The smaller 
(larger) the value of A the more the incipient donor stars are 
weighted toward higher (lower) mass. The contributions to 
the ULX population come from, in decreasing order of im- 
portance, black-hole binaries with: (i) high-mass compan- 
ions, (ii) mid- to high-mass donors on the giant branch, and 
(iii) lower-mass companions that lead to transient behavior. 
Second, our results are sensitive to the factor by which we 
allow for super-Eddington luminosities. In most of our cal- 
culations we have used a factor of 10 times the usual values 
of I/Edd (as estimated by Begelman 2002) . This substantially 
helps in the production of ULXs. However, we find that any 
further increase in this enhancement factor does not signifi- 
cantly increase the number of ULXs since the systems then 
become limited by the available mass-transfer rates. 

We find generally encouraging agreement between the 
ULX populations that we are able to generate and the 
observations — at least for the ULX luminosity function be- 
low ~ 2 X lO^'' ergs s~^. Some of our models, in particular 
Model B, yield substantial numbers of ULXs for times up to 
30 Myr after a star formation event, and a few such sources 
for times up to 100 Myr. We have found, however, by com- 
paring our production rate of ULXs and the concomitant 
values of ix.tot with the empirical relation between the star 
formation rates and I/x.tot derived by Gilfanov et al. (2004) 
and Ranalli et al. (2003) that our adopted black-hole bi- 
nary formation rate of 10~® yr~^ for a core-collapse SNe 
rate of 0.01 yr~^ is likely too high by a factor of ~ 3. Our 
calculated evolution of the luminosity function after a star 
formation event can explain spatial offsets between the loca- 
tion of ULXs and current-epoch star forming regions (e.g., 
Zezas & Fabbiano 2002; Fabbiano & White 2004, and ref- 
erences therein), provided that the wave of star formation 
passed from the former to the latter. This is a possible alter- 
native explanation to the one which requires that the ULXs 
are given natal kicks and ejected from the star clusters in 
which they are born (e.g., Zezas & Fabbiano 2002). 

Throughout this work all of the cited luminosities in the 
figures and table have been bolometric. On the other hand, 
most of the literature on ULXs cites Lx in the 2 — 10 keV 
band due to the sensitivity range of X-ray telescopes. It is 
difficult to estimate what fraction, /, of the bolometric X-ray 
luminosity emerges in the 2 — 10 keV band precisely because 
very few measurements extend to substantially lower or 
higher energies. Nothwithstanding this limitation, we have 
used various spectral shapes for ULXs reported in the liter- 
ature to estimate /. We find values of / between 0.1 and 0.5 
with a most probable value of ~0.3. This range for / should 
hold unless a large fraction of the luminosity is emitted in a 
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low-temperature multi-component disk spectrum that falls 
largely outside the 2—10 keV range. (For a related discussion 
see Appendix A2 of Portegies Zwart, Dewi, & Maccarone 
2004.) We can foresee two distinct possibilities for the fac- 
tor /: (1) it falls within a fairly narrow range (e.g., near 
~0.3) for most ULXs, or (2) it varies from source to source 
or temporally for individual sources. In the first instance, 
the entire luminosity functions that we generate could be 
shifted downward in luminosity such that, e.g., the peak of 
the luminosity function for model B drops from ~ 2 x 10^^ 
ergs s~^ to ~ 6 X 10'^** ergs s~^, and the highest (and rarest) 
luminosity sources of ~ 2 X 10''° ergs s'^ fall to ~ 6 x 10='*' 
ergs s"'^. By contrast, in the second case, if / varies from low 
to high values (among sources or temporally) then our lumi- 
nosity functions might be not be so seriously shifted to lower 
luminosities, but would yield fewer luminous sources over- 
all. This, in turn, would actually help the "overproduction" 
problem discussed in §3. If the net result of smaller values of 
/ is to shift our simulated luminosity functions to substan- 
tially lower luminosities, then this suggests two possibilities: 
(1) the basic results we have generated are correct, but then 
the "Begelman factor" would have to be as high as ~30 in- 
stead of 10 for black- hole accretors with masses of ~ 10 Mq, 
or (2) the most luminous ULXs with L^(2 - 10 keV) > 10"*° 
ergs s~^ indeed represent a different class of objects, e.g., 
IMBHs. More work is needed both observationally, to deter- 
mine / empirically with greater certainty, and theoretically, 
especially regarding the photon bubble instability model to 
ascertain if super-Eddington factors of ~ 30 are possible. 

In this work we have not included the contributions 
from very luminous neutron star binaries (e.g., HMXBs re- 
sembling SMC X-1 and LMC X-4). Such systems are Ukely 
to be very short lived (e.g., ~ 10^ yr; see e.g., Levine et al. 
1991; 1993; 2000) due to the tidal interactions that drive the 
Roche lobe rapidly through the atmosphere of the primary. 
However, such neutron star binaries are much easier to form 
than black-hole binaries, and therefore their contributions 
to the luminosity function between ~ lO^'^ and 5 x 10^^ ergs 
s~^ should be substantial. In any study attempting to com- 
pute the entire luminosity function, the luminous neutron 
star systems clearly need to be included (e.g., Belczynski et 
al. 2004). The intention of the present study, however, was 
limited to evaluating the contribution of stellar-mass black- 
hole binaries to the ULX population. 

One further caveat is in order with regard to the suc- 
cess of our model for the ULXs. Most of the ULXs that we 
generate in our population synthesis have high-mass donors 
and are not transients. By contrast, many of the black-hole 
binaries discovered and studied in our Galaxy are transients 
whose time-averaged mass-transfer rate is well below Ed- 
dington and whose companions are rather low in mass (i.e., 
< 1 Mq) (see, e.g., McClintock & Remillard 2004). This ap- 
parent discrepancy between the two populations was dis- 
cussed in PRH, and hinges on the fact that it is difficult, 
in the context of the present model, to produce black hole 
binaries with low initial donor masses. This problem arises 
due to the fact that a secondary of < 2 Mq tends to merge 
with the core of the massive black-hole progenitor instead 
of ejecting the common envelope. Also, as demonstrated in 
PRH, black-hole binaries with initial donor masses > 2 Mq 
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Figure 11. Potential X-ray luminosities as a function of time 
since the beginning of mass transfer (offset by 10^ yr for clar- 
ity) for 4 black-hole X-ray binary evolution sequences. Each color 
corresponds to a different mass or evolutionary state of the donor 
star when mass transfer commences (as indicated). In all cases 
the initial mass of the black hole is 5 Mq. The black and red 
curves are for donor stars with initial masses of 6 Mq , while the 
green and blue curves are for donor stars with initial masses of 
9 Mq . For the black and green curves the donor star commences 
mass transfer in the middle of the main sequence (case A mass 
transfer). The red and blue curves correspond to the case where 
the donor star has evolved off the main sequence and traverses the 
Hertzsprung gap at the start of mass transfer (early case B mass 
transfer). In all four calculations, thermal timescale phases last 
of order 10^ yr. Since these involve somewhat lower-mass donor 
stars than the main contributors to our ULX candidates, they 
may be more common, which can compensate somewhat for the 
shorter duration of the high mass transfer rate episodes. 



can evolve into the type of compact systems often seen in 
the Galaxy — but only if there is a significant source of sys- 
temic orbital angular momentum loss, e.g., magnetic brak- 
ing, which has been assumed not to operate in stars with 
radiative outer envelopes (Parker 1955; Pylyser & Savonije 
1988). We note, however, that there is a subset of stars with 
radiative envelopes, the Ap/Bp stars, with masses up to 
~ 3 Mq that have strong magnetic fields (~ 10 kG) and 
long rotation periods (e.g., Hubrig et al. 2000). While the 
origin of the large magnetic fields in these stars and their re- 
lation to the rotation periods are presently not understood, 
their very existence may suggest that there could be a sub- 
set of stars with radiative envelopes where magnetic braking 
is operative. This is clearly a problem that warrants further 
consideration and modeling. 

King et al. (2001) discuss the possible importance of 
thermal timescale mass transfer (see, e.g., PRH) onto black 
holes as the driver of high rates of mass transfer in ULXs. We 
note, however, that most ULXs in our simulations are bina- 
ries where mass transfer is driven by the nuclear evolution of 
the massive donor, most commonly hydrogen burning in the 
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core, but in a fraction of systems ( < 5 %) by hydrogen sheU 
burning as the donor star ascends the giant branch. Even 
higher mass-transfer rates can indeed be attained during 
phases where the secondary is out of thermal equilibrium, 
either during thermal timescale mass transfer or when the 
secondary moves across the Hertzsprung gap (so-called early 
case B mass transfer). However, since the lifetimes of these 
phases are typically several orders of magnitude shorter than 
the phases where mass transfer is driven by the nuclear evo- 
lution, few ULXs are expected to be found in this phase 
(see the dashed curve in Fig. |^^. To illustrate this further 
we have done some additional binary calculations, shown in 
Fig. 111! for secondaries with initial masses of 6 and 9 Mq, 
respectively, in two evolutionary phases: in the middle of 
the main sequence (case A) and just after the secondary 
has left the main sequence and started to move across the 
Hertzsprung gap (early case B). In all calculations, the ini- 
tial mass of the black hole was taken to be 5 Mq. In the 
6 Mq case A sequence, the secondary always remains in 
thermal equilibrium and mass-transfer is entirely driven by 
the nuclear evolution of the core. For about ~ 10^ yr its 
potential X-ray luminosity (just) exceeds 10'^^ ergs s^^ and 
may be classified as a ULX during this phase. In contrast, 
the 9 Mq case A sequence experiences a phase of rapid, 
thermal timescale transfer, lasting ~ 2 x 10^ yr. During this 
phase the potential X-ray luminosity reaches a peak exceed- 
ing 10''^ ergs s~^. After the thermal timescale phase, mass 
transfer continues but now on a nuclear timescale for the 
next 3 X lO'^ yr at a rate that is not sufficient to power a ULX. 
Both case A systems appear as ULXs again as they ascend 
the giant branch where their evolution is driven by hydrogen 
shell burning. In the case B sequences, where mass transfer is 
driven by the thermal timescale evolution of the secondary 
across the Hertzsprung gap, the mass-transfer phases last 
only ~ 10^ yr, but reach peak mass-transfer rates of 10~* 
and 8x 10"* Mq yr~^ for the 6 and 9 Mq sequences, respec- 
tively, and corresponding peak potential X-ray luminosities 
of order and above 10*^ ergs s~^. Note that in these latter 
sequences, the mass-transfer rates, and hence potential X- 
ray luminosities, drop by 1.5 orders of magnitude after the 
mass-ratio has been reversed and the orbit starts to expand 
rapidly. During this latter phase, the potential X-ray lumi- 
nosities remain in the ULX regime, but decrease gradually 
(the moderately sloping portions of the curves in Fig llH . 

These sequences fFig. lllH demonstrate that during ther- 
mal timescale mass-transfer phases, very high mass-transfer 
rates are attained, providing potentially enough fuel to 
power even the most luminous ULXs, but that the dura- 
tion of these phases, which obviously scales as the inverse of 
the characteristic mass transfer rate, is short compared to 
the duration of phases where mass transfer is driven by the 
nuclear evolution of the secondary. A compensating factor is 
that lower-mass secondaries can in principle produce mass- 
transfer rates in the ULX regime during thermal timescale 
phases (onto low-mass black holes) that they could not oth- 

^ This is generally true for mass transfer from a more massive 
star to a less massive star and, as is well known, explains, e.g., 
the paucity of (3 Lyrae systems compared to Algol-type systems, 
believed to be their descendants (for a general discussion, see 
Paczynksi 1971). 



erwise do during the nuclear-driven phases. This effect may 
somewhat, but not drastically increase the relative number 
of thermal timescale systems compared to nuclear timescale 
systems. 

Transient source behavior, due to the thermal- 
ionization disk instability (e.g., van Paradijs 1996; King, 
Kolb, & Burderi 1996; Dubus, Hameury, & Lasota 2001) 
has been invoked as a way of producing higher mass trans- 
fer rates, albeit only for concomitantly reduced fractional 
times of the total binary mass transfer phase. We have found 
that when we include transient source behavior, which oc- 
curs predominantly in systems with lower average luminosity 
(e.g., < lO^'^ ergs s~^), there is essentially no change in the 
luminosity function for times < 100 Myr; thereafter, there 
is a modest enhancement of high-La, sources due to tran- 
sient source behavior, but the bulk of these objects have 
luminosities ^5 x 10^* ergs s~^ (see Fig.|l] panel [d]). 

The giant phase of the donor star is also a natural 
mechanism to consider for driving mass transfer rates that 
could enhance the ULX population. The mass transfer phase 
driven by the ascent of the donor up the giant branch is dra- 
matically evident in Fig.Qas "spikes" in near the end of 
the evolution. Indeed the mass transfer rate increases during 
this phase by 1 - 2 orders of magnitude. However, this phase 
is relatively short lived and tends not to contribute dramati- 
cally to the population of ULXs. Fig. 0|shows the integrated 
contribution of the giant branch phase of mass transfer to 
the ULX population — these are seen as a red ridge lying 
~1.5 orders of magnitude above the peak in the luminosity 
function for times > 15 Myr. One could argue that in galax- 
ies where the Chandra sensitivity limit is ~ 10''^ ergs s~^, 
these (giant branch) sources would be the only ones visi- 
ble at tcv ^ 15 Myr. This is true, but then the issue arises 
as to what absolute numbers of sources are predicted. As 
can be seen from Fig. |S]the contributions from giant-branch 
donor sources do show up as a modestly significant feature 
for > 10**" ergs s"^ at ~30 Myr. 

Our models also have some predictive power concerning 
the optical appearance of the ULX black-hole binaries. We 
have computed tracks for our binary evolution models in the 
HR diagram, taking into account both the contribution from 
the donor star and the accretion disk. Light from the disk, 
in turn, results from reprocessing of X-radiation and from 
viscous heating. We computed the effective temperature of 
the disk as a function of radius, r, from the following simple 
expression that we derived: 

where x is the radial distance in units of r^in, the inner ra- 
dius of the disk, a is the X-ray albedo of the disk (which we 
take to be ~0.7), and ^' is defined by ^' = ^(rmin/fmax)^^^, 
where, in turn, ^ is the half thickness of the disk, h, at rmax, 
in units of rmax. In this formulation h{r) — £,r'^^'^ra^J^ , and 
we take ^ = 0.1 as an illustrative value corresponding to a 
full angular thickness of the disk equal to ~ 12°. For the 
disk-heating problem only, we have taken the inner edge 
of the accretion disk to be at QGM^yi/ <? regardless of the 
spin of the black hole, and we have neglected the factor 
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Figure 12. Evolution tracks in the HR diagram are shown for 4 
of our models; these are for donors stars that had initial masses 
of 4, 6, 8, and 10 Mq and were on the zero-age main sequence at 
the start of mass transfer. Each binary system is represented by 
two different evolution tracks; the thin curve is the contribution 
from the donor star alone, while the thicker curve (of the same 
color) describes the track for the total system light - both disk 
and donor star. 

^1 — \Awn7^^ found in the Shakura & Sunyaev (1973) 

solutions. This latter approximation afltects only the a::"'^ 
term which, in any case, contributes little to the optical flux 
from the disk. The contributions to the B and V bands from 
the donor stars were estimated from their bolometric mag- 
nitude and effective temperature using conversion factors 
taken from Reed (1998). 

Illustrative results for 4 evolution tracks in the HR di- 
agram are shown in Fig. 1121 these are for donors stars that 
had initial masses of 4, 6, 8, and 10 Mq and were on the zero- 
age main sequence at the start of mass transfer. Each binary 
system is represented by two different evolution tracks in 
the HR diagram; the thin curve is the contribution from the 
donor star alone, while the thicker curve (of the same color) 
describes the track for the total system light - both disk and 
donor star. In each case, it is clear that during the early part 
of the binary evolution where the donor star is still on the 
main sequence, the contribution from the heated accretion 
disk leads to a brighter optical counterpart in the V band 
by about 1 magnitude. However, during the short-lived gi- 
ant branch phase, when the X-ray luminosity is very much 
larger, the V-band magnitudes are enhanced by ~ 4 mag- 
nitudes and the system colors are considerably bluer. Thus, 
we conclude that there is a substantial contribution to both 
the color and magnitude of the optical counterparts of the 
black-hole ULXs from the accretion disk. Moreover, we note 
that the contribution of the donor star to the optical light 
of the system may be considerably less than might be ex- 
pected from a star of the same initial mass, especially if the 
mass transfer commences early in the evolution of the donor 
star. These facts should be taken into account by observers 
trying to characterize the "donor star" from its location in 
an HR diagram (M. PakuU, private communication). 



Finally, we briefly discuss the issue of ULXs in ellipti- 
cal galaxies (e.g., Angelini, Loewenstein, & Mushotzky 2001; 
Colbert & Ptak 2002; Jeltema et al. 2003; Ptak & Colbert 
2004). The presence of true ULXs with certified X-ray lu- 
minosities exceeding 10'^^ ergs s~^ would be very difficult 
to explain within the context our black hole binary model 
unless the elliptical galaxy had a recent (e.g., ^ 3 x 10* yr) 
merger and star formation cycle. Recently Irwin, Bregman, 
& Athey (2004) have shown, from a sample of 28 elliptical 
and SO galaxies observed with Chandra, that the number of 
sources with L^, >2 x 10^^ ergs s~^ was equal to the num- 
ber expected from background or foreground objects. Thus, 
they conclude that, with only the exception of two ULXs 
in globular clusters within NGC 1399, ULXs "are gener- 
ally not found within old stellar systems". Nonetheless, the 
ULXs found in early-type galaxies with 10^® < < 2 x 10^^ 
ergs s~^ are not naturally explained within the context of 
our model. 

In summary, our calculations indicate that with a plau- 
sible set of assumptions, a majority of the ULXs with 
Lx ^ lO''" ergs s~^ in spiral galaxies can be understood 
in terms of binary systems containing stellar-mass black 
holes. The systems can be evolved theoretically from the 
primordial binary phase through the ULX phase — and the 
absolute numbers of such systems in the populations that 
are computed are in agreement with the observations to 
order-of-magnitude accuracy. By contrast, the same claims 
cannot yet be made for the model of ULXs that invokes 
intermediate-mass black holes. Such a model, of course, has 
the advantages that (1) the Eddington limit presents no dif- 
ficulties, and (2) the spectra of some of the ULXs which 
exhibit low inner-disk temperatures (see, e.g., Miller et al. 
2003; Cropper et al. 2004) may be a natural consequence 
of more massive black holes. However, in the context of the 
IMBH model, ideas as to how such objects are formed, ac- 
quire binary companions, and are "fed" at the requisite ac- 
cretion rates are only just beginning to be explored (see, 
e.g., Portegies Zwart et al. 2004; Portegies Zwart, Dewi, & 
Maccarone 2004; and references therein). 
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